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Abstract. A new method for the determination of electric signal time-shifts is introduced. As the 
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Kolmogorov-Smirnov test, it is based on the comparison of the cumulative distribution functions of the 
reference signal with the test signal. This method is very fast and thus well suited for on-line applications. 
It is robust to noise and its performances in terms of precision are excellent for shift times ranging from a 
fraction to several sample durations. 



> 

PACS. 29.40. Gx Tracking and position-sensitive detectors - 29.30.Kv X- and gamma-ray spectroscopy - 



(N 



07.50.Qx Signal processing electronics 



1 Introduction the characteristics of the particles. The properties of the 

o 
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particle actually detected are supposed to be the same 
Modern detection systems depend more and more on em- ^s those of the best-matched signal of the basis. However, 



^ bedded intelligence that allows the on-line determination random time jitters often alter detector signals, which can 
of the characteristics of the impinging particles [ll|2j[3l|4] . severely affect their comparison to basis signals even when 
In many cases, the full shape of the delivered signals is the time-shift is smaller than the sample duration |5]. 
first digitized and then compared to a basis of reference More generally, signal absolute-time alignment with re- 
signals. The basis signals reflect all the possible values of gpect to a reference signal or relative-time alignment (e.g. 
Correspondence to: Pierre.Desesquelles@in2p3.fr with respect to the signal set average) are important is- 
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sues in detector signal processing. In most cases, on-line In Section 2, we summarize the KS method and show 

time-correction techniques, such as constant fraction dis- why it is more sensitive to curve shifts than the usual chi- 

crimination, rely on the crossing of a threshold. Due to square method. An illustration of the method, using High 

the signal noise, and in many applications to the digitiza- Purity Germanium (HPGe) detector signals, is given in 

tion, these techniques do not allow time correction shorter Section 3. The conclusions are enumerated in Section 4. 
than the sample duration (intra-sample correction). The 



2 The Kolmogorov-Smirnov test 



poor-man solution to time jitters consists in a chi-square 
comparison of the detected signal with different versions of 

time-shifted reference signals. This solution is often inap- "^-^ Description of the method 

plicable as it may severely slow down the acquisition rate. 
In the following, we propose an algebraic determination 
of the time-shift based on the Kolmogorov-Smirnov (KS) 
distance. This method consists in measuring the time dis- 
tance between the detected signal and a reference signal. It 
is well adapted to on-line applications since it is very fast 
and it entails only a small memory space. It will be shown 
that the resolution of the KS method is excellent, with 
respect to the usual time alignment methods, for noisy 
signals and when the time-shift is a few sample durations 
or shorter. Therefore, for longer time shifts, it can be used 
to refine a first evaluation obtained by another method. 



The KS distance is, as the chi-square, a criterion to qualify 
the agreement between two distributions. Its implemen- 
tation is illustrated in Figs. [T] and [2] In the first figure 
a "theoretical" distribution is compared to "experimen- 
tal" points affected by statistical fluctuations. The KS 
distance between the two distributions is defined as the 
maximum distance I?ks between their cumulative distri- 
bution functions (see Fig. [2]). The cumulative distribution 
function of a probability density function / is defined as 
F{x) — f{x') dx' . We use the same term for discrete 
functions: F(j) = Yfi=i fii ^'^^^ fi = "i/^tot, where rii is 
the number of counts in the i"^ sample and A^tot the total 
The method is simple to implement and should be well number of counts. Cumulative distribution functions are 
suited to a large number of appHcations. However, it is not monotonically increasing from zero to one. 
the case for the measurement of gas chambers drift-times. The so-called chi-square distribution gives the proba- 
which could have been an important application. Indeed, bility that the measured un-normalized chi-square (Xm) 
the KS method requires well defined reference signals but superior to a given value. Such a distribution also exists for 
also that the detector signals be affected by uncorrelated the KS distance (Fig.jsj). In our case, if the differences be- 
white noise only. Gas chamber signals are altered by strong tween the theoretical curve of Fig. [T] and the experimental 
auto-correlated fluctuations due to the diffusion phenom- data are only due to the fluctuations, then the chi-square 
ena in the gas. has a 5% probability of being larger than 62 and the KS 
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Fig. 1. Comparison of a generic theoretical distribution (full 
line) with experimental count numbers (black dots) as a func- 
tion of the sample number (this figure shows an example made 
by hand only to illustrate the method, it does not present real 
data) . 
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Fig. 2. Comparison of the cumulative functions of the distri- 
butions from Fig. [l] The KS distance is the absolute value of 
the maximum difference between the two. 



distance has a 5% probability of being larger than 0.066. 
If the measured values are greater than about 70 and 0.08 
respectively then the theoretical curve is most probably 
invalidated by the experimental results. 



Fig. 3. Probability that the measured chi-square (left) and 
measured KS distance (right) are greater than a given value. 
The dotted lines correspond to a probability of 5%. The shape 
of the chi-square distribution depend only on the number of 
samples, while the shape of the KS distance distribution de- 
pends only on the total number of counts. 

The KS distance and the chi-square are not sensitive 
to the same types of discrepancies between the distribu- 
tions under comparison. It is visible from Fig. [T]that the 
theoretical curve is shifted to the left with respect to the 
experimental points. However, the curve crosses most of 
the error bars. Hence, the measured chi-square has a low 
value (Xm — 40.4) from which one could conclude that 
both distributions are in agreement. Indeed, as shown in 
Fig. [3] the probability to obtain a chi-square equal or 
greater than 40.4 is 67%. However, the KS distance for 
the same distributions is 0.084, which corresponds to a 
very low probability of 0.6 %. This shows that the theo- 
retical curve is most probably invalidated by the data. The 
KS distance is a test particularly sensitive to the shifts be- 
tween distributions. In the following, this well established 
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result of probability theory will be applied to determine where At is the time-shift between the test and the refer- 
the time-shifts of detector signals affected by noise and a ence cumulative distributions, 
random time jitter. 



2.2 Application to tlie determination of random 
time-shifts 

The cumulative function S is now obtained by integrating 
the signal s and normalizing the maximum to one: 



Sit) 



(1) 



El-sit') 

Our goal is not to measure precisely the KS distance be- 
tween the shifted test signal s and the reference signal sq, 
but to measure the time-shift in a fast way. The verti- 
cal distance between slightly shifted distributions is large 
where the slope is large, as can be seen for example in 
Fig. [2] since, to the first order, AS = (dS'/dt) At. The 
maximum slope of the cumulative distribution corresponds 
to the maximum of the signal (in the example of Fig. [l] 
^Om = 13). Therefore, the distance between the reference 
signal cumulative distribution Sq and the test signal cu- 
mulative distribution S will be evaluated at t = iom- 

Noting toc the central time of the sample which con- 
tains the reference signal maximum (see Fig. |4]), and tg 
the sample duration: 



(2) 



we have: 
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Fig. 4. Scheme of the determination of the time-shift using 
the Kolmogorov-Smirnov method. The full line curve is 5*0, 
the reference signal cumulative function (1 ns time bins) and 
the dotted curve is S, the test signal cumulative function (10 ns 
time bins). The bold curve is the third order Taylor expansion 
of the reference signal around tom. The estimation of the time- 
shift is the distance between the dotted and the bold curves at 



The value of the reference cumulative distribution func- 
tion can be evaluated using its Taylor polynomial [B]: 



Soitoc + At) — S'o(iom) 



(4) 



+ (<0c + At- tom) 



dSo 
dt 



+ ^ (ioc + At- tom)^ 



d^So 



dt^ 



+ ^ itoc + At- tom)^ ^ 



Sitoc) = Soitoc + At) , 



(3) 
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The cumulative distribution function is obtained by inte- statistical fluctuation on the cumulative distribution to 
grating the signal ([T]) thus, its first derivative is simply the the time-shift evaluation is minimized, 
signal itself: (dS'o/dt)(tom) = so(iom)- The second deriva- The complete third order solution reads: 
tive is equal to zero as it corresponds to the slope of the 

signal when it reaches its maximum. The third derivative (^i _ iy^) (1 + i\/3) ct 

^tKS — ;^ 1- ^Om — ^Oc , (8) 



can be estimated by: 



so(iom + 2) - 2so(^Om) + so(iom - 2) 

'^Om = , (5) 



where tg is the sample duration (1 ns in our case). Using 
Eqs. (|3|, ^ and ([s]), and simplifying the notations, one 
obtains: 



2r-on 

ith: 



a = [^J rg^(8 sg„j -I- 9 ro^{Sc - S'om)^) 

+ Zrl^{S,-S^^)Y . (9) 



3 Results 

3.1 The set of test signals 

= 5*0111 (6) 

-f (ioc + ^^KS ~ ^Om) som Order to illustrate the method, we apply it to large sets 

, 1 /, , /i, J. \3 of simulated signals with different amounts of noise and 
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In this equation, all the constants ^om, ^Oc S'om, som, J'om 
depend only on the reference signal and thus, can be pre- 
calculated off-line. Only these constants have to be saved 
in the computer memory, instead of the full signals. The 
value of 5c is calculated by integrating the test signal. 
Hence, Eq. ([6]) is a completely defined third order equation 
in Z\iKS • The first order solution (which is also the second 
order solution since (d^5o/di^)(tom) = 0) reads: 



different time-shifts. These signals were generated using 
the MGS code [ZIIHIIS], that simulates HPGe signals [TU]. 
A large variety of rise times and pulse shapes is obtained. 
The sample duration of the simulated signal is 1 ns. As the 
experimental signals, they are grouped into 10 ns samples. 
Eventually, white noise is added to the signals. The refer- 
ence signals are the 1 ns plane signals and the test signals 
are the 10 ns noisy and shifted signals (respectively in full 
and dotted lines in the example of Fig.|5|. The KS method 
applies to positive distributions, thus, in the case of non- 

Zi^KS = ^'^ ~ "^""^ + ^0 ~ ^0 (7) 

som ™ positive electric pulses, their absolute values must be con- 

The derivative of the resulting time-shift with respect to sidered. In the specific case of segmented HPGe detectors, 

the test signal is (dAtyjil dS^ = l/'^Om- As som is the there are two types of signals: transient signals, and sig- 

maximum of the signal, the error propagation from the nals from the hit segments which are charge-integrated, 
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SO that their shape is similar to a cumulative distribution 
function. In the latter case, the KS method can be applied 
directly to the signal, after normalization of its amplitude 
to 1. Unfortunately, in this case, one does not benefit from 
the smoothing property of signal integration. 




100 200 300 400 100 200 300 400 

f (ns) t (ns) 



Fig. 5. One example of signal (left column) and the corre- 
sponding cumulative distribution function (right column). The 
reference signal (1 ns samples) is represented in full line and the 
test signals (10 ns samples) are represented in dotted line. The 
test signals are shifted by 10 ns with respect to the reference 
signal. The amount of noise is indicated for each line. 



In the following, we do not consider the problem of 
finding the basis reference signal that corresponds to the 
test signal. This problem is crucial in many applications. 
For discussion of the methods to determine the right ref- 
erence signal, in the case of HPGe detectors, see [51|5J[ni 
[TOjITT] and references therein. 



3.2 KS method 

The application of the first order and third order cor- 
rections to noisy and shifted test signals is illustrated in 
Fig. [6] The plots show the error on the correction as a 
function of the imposed time-shift for different values of 
the signal-to-noise ratio. The down-bendings of the curves 
are due to the truncation of the third and fourth order 
terms in the Taylor expansion. 

Inside the first 10 ns sample, both corrections are un- 
biased whatever the noise amount. The third order cor- 
rection is still unbiased at the end of the second bin. Be- 
yond the second bin, a small negative drift is visible (the 
time-shift is slightly underestimated) . The error bars show 
that the method reaches its limit {Atyis — /ii ~ 10 ns) at 
SNR — 3 dB which is a usual minimum value of the 
signal-to-noise ratio when extracting information from a 
noisy signal. 

3.3 Comparison with other methods 

A large number of methods for the determination of time- 
shifts are available. Electronic or software CFD are 
fast but very sensitive to noise and, except for very high 
signal-to-noise ratios, not adapted to intra-sample correc- 
tion. The chi-square or residue methods consist in the ex- 
tensive comparison of the test signal with different, reg- 
ularly time-shifted, versions of the reference signal jTT] . 
This procedures allow a comparable time resolution but 
are about two orders of magnitude slower than the solu- 
tion we propose and entails the handling of bases of shifted 
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Fig. 6. Error on the time-shift determination using the 
Kolmogorov-Smirnov method. The black dots correspond to 
the first order correction and the open dots correspond to the 
third order correction. The error bars give the one-sigma width 
(one thousand test signals per dot). 

reference signals that may be very large. Another popular 
technic consists in calculating the centroid time [13 of the 
detected pulse : iccn = J2t^ ^i^)- ^^'^ computa- 

tion time is the same as for the KS method (Eqs. ^ and 
0). The comparison of the absolute value of the error 
fot both methods is presented in Fig. [Tj The KS method 



appears to be less sensitive to the noise and its resolution 
is much better for time-shifts smaller than a few sample 
durations. 
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Fig. 7. Comparison of the average absolute value of the error 
on the time-shift determination for the first order KS method 
(black dots) and the method of the centroid (open dots) as a 
function of the signal-to-noise ratio. The simulated time-shift 
ranges from 10 ns to 80 ns. The results obtained with the 
centroid method are almost independent of the imposed time- 
shift. The sample duration is 10 ns 



4 Conclusion 

We have presented a method, based on the Kolmogorov- 
Smirnov test, to correct signals corrupted by time jitters. 
This method appears to be very efficient for the precise 
correction of time-shifts up to a few bin times. For SNR 
larger than 20 dB, the bias for intra-sample correction is 
less than one tenth of the sample duration. The good pre- 
cision of the corrections is connected to two features of 
the method. On the one hand, the time-shift is calculated 
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from the cumulative distribution function of the signal (so 
that the noise fluctuations are lowered by the integration, 
see Fig. [s]), and, on the other hand, it is evaluated where 
the pulse amplitude is maximum so that the error propa- 
gation is minimized. 

The first order correction is extremely fast since it en- 
tails only the integration of the signal and the calculation 
of a first degree polynomial, Eq. ([t]) . The third order cor- 
rection is non- linear complex, Eq. Its bias is smaller 
than in the first order case, but, as it includes the calcu- 
lation of higher powers, the uncertainties are larger. 

In our AG ATA example, the reference signals are sim- 
ulated. The method can also be used for signal alignment. 
In this case, the reference signal could simply be the mean 
of the measured signal set. 

Finally, due to its simplicity, the method is versatile 
whereas it is probably not suited to any type of signals 
(for example when the signals cannot be corrected from 
baseline / DC offsets). In the case of HPGe detector sig- 
nals, its reliability is comparable with other methods for 
long time-shifts, while, in the range of one tenth to a few 
sample durations the results, in terms of resolution, ro- 
bustness to the noise and computing time, are excellent. 
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